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Abstract. We have discovered with the Wide Field Cam- 
eras on board BeppoSAX the weak transient X-ray source 
SAX J2239. 3+6116 whose position coincides with that of 
4U 2238+60/3 A 2237+608 and is close to that of the 
fast transient AT 2238+60 and the unidentified EGRET 
source 3EG 2227+6122. The data suggest that the source 
exhibits outbursts that last for a few weeks and peak to a 
flux of 4 x 10~ 10 erg cm~ 2 s _1 (2-10 keV) at maximum. 
During the peak the X-ray spectrum is hard with a pho- 
ton index of —1.1 ± 0.1. Follow-up observations with the 
Narrow-Field Instruments on the same platform revealed 
a quiescent emission level that is 10 3 times less. Searches 
through the data archive of the All-Sky Monitor on RXTE 
result in the recognition of five outbursts in total from this 
source during 1996-1999, with a regular interval time of 
262 days. Optical observations with the KPNO 2.1 m tele- 
scope provide a likely optical counterpart. It is a B0 V to 
B2 III star with broadened emission lines at an approxi- 
mate distance of 4.4 kpc. The distance implies a 2-10 keV 
luminosity in the range from 10 33 to 10 36 erg s _1 . The 
evidence suggests that SAX J2239. 3+6116 is a Be X-ray 
binary with an orbital period of 262 days. 

Key words: stars: individual: SAX J2239.3+6116, 
4U 2238+60, M 2233+595, 3A 2237+608, AT 2238+584, 
3EG J2227+6122 



1. Introduction 

The Wide Field Cameras (WFCs) on board the BeppoSAX 
satellite carry out an extensive monitoring program of the 
whole X-ray sky. Various kinds of transient X-ray sources 
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are detected during these observations on time scales be- 
tween a fraction of a second and months. We here discuss 
one case which we have followed up with the more sensitive 
Narrow-Field Instruments (NFI) on board BeppoSAX and 
with the Kitt Peak National Observatory (KPNO) 2.1 m 
optical telescope, and whose signature we searched for in 
data from the All-Sky Monitor (ASM) on the Rossi X-ray 
Timing Explorer (RXTE) and in occupation data from 
the Burst And Transient Source Experiment (BATSE) on 
the Compton Gamma-Ray Observatory. In all observa- 
tions we have interesting detections. We describe analyses 
of WFC (Sect. |), NFI (Sect, g), ASM (Sect. |), BATSE 
(Sect. ||) and optical data (Sect. ^|), point to possible asso- 
ciations of previously found sources of high-energy emis- 
sion in the field (Sect. [?]), and discuss the possible nature 
of the source (Sect. ||). 

2. WFC observations 

The WFCs (Jager et al. 1997) on the BeppoSAX satel- 
lite (Boella et al. 1997a) carry out a monitoring program 
where they point at arbitrary positions in the sky. This 
program is the 'secondary mode' program whereby the 
WFC pointings are dictated by those of the Narrow Field 
Instruments (NFI) on the same satellite and the con- 
straints of the satellite attitude with respect to the direc- 
tion of the sun. The secondary mode observations entail 
roughly 90% of all observations. 

There are two identical WFCs. They each have 40 x 
40 square degrees field of views (full width to zero re- 
sponse) with ~ 5' angular resolution. The bandpass is 2 
to 26 keV. The WFCs have been active continuously since 
mid-1996, except for a 1-month period in May 1998 when 
unit 1 was turned off, and a 3 week period in early 1999 
when both cameras were turned off. During mid-May to 
mid-August 1997, the satellite was in a standby mode, 
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with both cameras turned on but without sufficient stabi- 
lization for sensitive imaging. 

During 103.4 d on-source time with a net exposure 
time of 27.8 d, the secondary observations covered the 
field around the supernova remnant Cas A. This exposure 
is roughly uniformly distributed in time. The latest data 
that we report here were obtained on Dec. 15, 1999. Ex- 
cept for Cas A, there are no bright X-ray sources present 
in this field. Compared to the background levels in the 
WFCs, Cas A is not so bright. The sensitivity is, there- 
fore, close to optimum in the field. During two instances, 
a relatively faint transient turned up in the field, at an 
angular distance of 6.0° from Cas A. They were on March 
4, 1997, and May 8, 1999. We note that sometimes Cyg X- 
2 is in the same field of view as this transient, but at an 
angular distance of 24.5° which implies it does not de- 
grade the signal of the transient. The celestial position 
of the WFC-detected transient is a 2 ooo.o = 22 h 39 m 22.8 s , 
^2000.0 = +61°16'47" with a 99% confidence error radius 
of 1.'8. 

The detection in March, 1997, was the strongest WFC 
detection obtained so far. In Fig. we present a light 
curve with a resolution of about a day. The period over 
which we detected activity is about 1 week. The variation 
in brightness during this event seems to suggest that it 
was not active for much longer. 

A light curve of the second WFC-detected outburst 
is presented in Fig. ^. The measured peak flux of this 
outburst is about a factor of 3 lower than that of the first 
outburst but this could very well be a sampling effect. The 
source was above the detection threshold for one week but 
due to an unfavorable sampling it could just as well have 
been active for 50 days. 

We have extracted the spectrum from the observation 
at the peak of the first outburst and modeled it with a 
simple power law plus absorption due to interstellar gas 
of cosmic abundances (according to Morrison & McCam- 
mon 1983). We kept the value for the hydrogen column 
density A^h fixed at 1.0 x 10 22 cm" 2 , as found from HI 
maps (Dickey & Lockman 1990). The fit was satisfactory, 
with Xr — 1-22 (26 dof). The photon index is quite hard 
at -1.1 ± 0.1. The spectrum is shown in Fig. |. The 2-10 
keV flux is 3.3 x 10~ 10 erg cm~ 2 s _1 which is 0.016 times 
the flux of the Crab in the same bandpass. In the energy 
range 2 to 26 keV the flux is 1.0 x 10~ 9 erg cm~ 2 s _1 or 
0.033 times that of the Crab. 

In Fig. ||, we present a light curve with a relatively high 
time resolution when the flux was at maximum during 
the first outburst. The source appears to be variable but 
no coherent signal could be found in this data set. The 
90%-confidence upper limit on the amplitude of a periodic 
sinusoidal variation is 10% for periods between 1 and 10 3 s. 

We have analyzed WFC data taken on December 13- 
15, 1999. This is 2 weeks after the optical measurements 
described in Sect, pi No X-ray detection was found. The 
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Fig. 1. Light curve of the first WFC-detected outburst of 
SAX J2239. 3+6116 in March, 1997, for 2 to 26 keV with 
a time resolution set by the duration of the measurements 
(these are of order 1 d). The arrows indicate 3rj upper 
limits. These limits vary considerably because of varying 
exposure times and varying exposed detector areas. March 
6.877, 1997, is MJD 50513.877. 
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Fig. 2. Light curve of the second WFC-detected outburst 
of SAX J2239.3+6116 for 2 to 26 keV. May 8.241, 1999, 
is MJD 51306.241. 

upper limit on the X-ray flux is 6xl0 -11 erg cm -2 s _1 
(2-10 keV). 

3. NFI observations 

Observations were made with the NFI on Nov. 24.46- 
25.04 (UT), 1998, which is 627 d after the first and 
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Fig. 3. Top panel: WFC-measured spectrum of SAX 
J2239.3+6116 (crosses) and the fitted power law model 
(histogram). Bottom panel: residuals between data and 
model in units of a. 
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Fig. 4. 2 to 26 keV light curve of SAX J2239.3+6116 on 
March 6, 1997. The time resolution is 75 s. 



165 d before the second WFC detection. The NFI in- 
clude the low-energy concentrator/spectrometer (LECS; 
0.1-10 keV; Parmar et al. 1997), the medium-energy con- 
centrator/spectrometer (MECS; 1.8-10 keV; Boella et al. 
1997b), the high-pressure gas scintillator proportional 
counter (HP-GSPC; 4-120 keV; Manzo et al. 1997), and 



100 



2CC 



•300 



400 



Fig. 5. MECS image of observation on SAX 
J2239.3+6116. MECS units 2 and 3 data are com- 
bined. 



the phoswich detector system (PDS; 15-200 keV, Frontera 
et al. 1997). The first two instruments are imaging de- 
vices. The effective exposure times on SAX J2239. 3+6116 
are 11.9 ks (LECS), 20.1 ks (MECS) and 10.0 ks (PDS). 
The HP-GSPC was turned off for operational reasons. 

Only the MECS data show a detection of 
SAX J2239. 3+6116 in an otherwise empty field of 
view (see Fig. |J). The position of the source is 
a 20 oo.o = 22 h 39 m 18.5 s , <5 2000 .o - +61°16'20" with a 
90% confidence error radius of 1 '0. This is 0'.7 from 
the WFC centroid position and the two positions are 
fully consistent (see Fig. ||). The chance probability 
of a background point source with a flux in excess of 
10 -13 erg cm~ 2 s _1 in the WFC error region is 0.3% 
(based on the counts by Giommi et al. 1999). We reduced 
the spectrum of the point source from the MECS data 
through application of the maximum likelihood method 
(e.g., In 't Zand et al. 2000) and tried to model it with 
absorbed single-component models like a power-law 
function and thermal bremsstrahlung, while fixing Nn 
to lxlO 22 cm~ 2 . Only black-body radiation gave a 
satisfactory fit to the data of a single-component model, 
with Xi = 1-74 for 7 dof. The color temperature is 
1.2 ± 0.2 keV and the emission area at 4.4 kpc (see 
Sect. H) for isotropic emission is equal that of a sphere 
with a radius of 0.6 ± 0.4 km. The 2-10 keV flux is 
5 x 10~ 13 erg cm~ 2 s _1 which is ~ 10 3 times fainter than 
the peak flux of the March 1997 outburst. 
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Fig. 6. (a) ASM light curve at 4 d time resolution, the 
arrows point to the predicted times of outbursts for pe- 
riodicity of 262 d; (b) 9 statistic (Stellingwerf 1978) as a 
function of test period; (c) folded light curve at the best 
fit period. The non-zero flux outside the outburst points 
to an instrumental bias level of about 0.5 mCrab. 

4. ASM observations 

Given the NFI position, we have searched for signatures 
of SAX J2239.3+6116 in archival ASM data. The ASM 
is fully operational since March 1996 and monitors each 
position on the sky in 2 to 12 keV during 90 sec snap- 
shots with a frequency of 5 to 10 times a day at a sensi- 
tivity of ~ 10 mCrab per day of observervations on un- 
crowded fields (Levine et al. 1996). The lightcurve for SAX 
J2239. 3+6116 at 4-day resolution is presented in Fig. ||a. 
There is the suggestion for a detection during 5 instances, 
at regular interval times of about 262 d. We tested this pe- 
riodicity by first filtering out the data within 10 d of clos- 
est aproaches to the Sun and all data after MJD 51450 
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Fig. 7. BATSE light curve of the outburst of SAX 
J2239. 3+6116 in March, 1997, for 20-100 keV with a one 
day time resolution. 
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Fig. 8. Top panel: BATSE-measured spectrum of SAX 
J2239.3+6116(crosses) and the fitted power-law model 
(histogram). Bottom panel: residuals between data and 
model in units of sigma. 



because there does not appear to be an outburst there, 
and then calculating the variance statistic 9 as defined by 
Stellingwerf (1978) for a range of test periods, see Fig. ^b. 
The resulting period is 262 ± 5 d. The epoch for peak flux 
is MJD 50786. The predicted times of outbursts are indi- 
cated in again Fig. ^a. Only the last predicted outburst 
does not appear to have materialized. The two WFC de- 
tections synchronize with the 2nd and 5th outburst. A 
folded light curve (Fig. ||c) shows an average outburst 
profile which lasts ~15 days. The average peak flux is 
4 mCrab. 



5. BATSE observations 

Due to the hardness of the spectrum measured with the 
WFC a search of the BATSE data was made to see if 
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SAX J2239. 3+6116 was detectable in the hard x-rays. The 
BATSE experiment onboard the Compton Gamma-Ray 
Observatory (Fishman et al. 1989) using the Large Area 
Detectors (LADs) can monitor the whole sky almost con- 
tinuously in the energy range of 20 keV to 2 MeV with 
a typical daily 3cr sensitivity of better than 100 mCrab. 
Detector counting rates with 16 energy channel energy res- 
olution and a timing resolution of 2.048 seconds (CONT 
data) are used for our data analysis. 

To produce the SAX J2239.3+6116 light curve, single 
step occultation data were taken using a standard Earth 
occultation analysis technique used for monitoring hard X- 
ray sources (Harmon et al. 1992). Interference from known 
bright sources was removed. The single occultation step 
data were then fit with a power law with a -2.0 photon 
index to determine daily flux measurements in the 20-100 
keV band. A time period of 30 days centered on each out- 
burst was examined for emission. Various time averages of 
the data were performed to search for a possible signal in 
the data. 

During the time of the March 1997 outburst we find 
evidence of a weak two day outburst in the BATSE data 
(see Fig. 0). The 20-100 flux for highest daily average cor- 
responds to 0.080 times that of the Crab in the same band- 
pass. It appears that there is a ~3 day delay between the 
peak in the BATSE data and that detected in the WFC. 
For any of the other outbursts predicted throughout the 
BATSE operational period we do not find evidence for 
hard x-ray emission, but we can not rule out the presence 
of a source with a flux 0.050 times that of the Crab or less 
being present in the data. 

We have extracted spectra from the BATSE data for 
the March 1997 outburst. During this time period the 
source was measured in three detectors. Spectra and re- 
sponses for each detector were created and a joint fit was 
made with a simple power law. The fit was satisfactory, 
with a reduced chi-squared = 0.90 (40 dof). The photon 
index is —2.2 ± 0.3. In Fig. || is plot the spectrum in which 
the three data sets have been grouped for plotting pur- 
poses. The 20-100 keV flux is 1.4 x 10~ 9 erg cm -2 s -1 
which is 0.080 times that of the Crab in the same band- 
pass. Very little flux is detected above 100 keV. 

6. Optical observations 

We looked for the optical counterpart of 
SAX J2239. 3+6116 by obtaining spectra of the three 
brightest stars in the NFI error circle on 1999 Dec. 2 and 
3 UT, using the KPNO 2.1 m telescope and Goldcam 
spectrograph. The date is 50 d before a predicted outburst 
(the last indicated in Fig. ||a). The spectra covered the 
wavelength range 3900-7500 A at 5 A resolution. The 
second brightest star within the error circle (Fig. |^) is 
the only one that has broad emission lines and other 
spectral features that are characteristic of a reddened 
Be star. Be stars are often hard X-ray sources because 
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Fig. 9. Finding chart for the optical counterpart of 
SAX J2239. 3+6116 from the Digitized Palomar Obser- 
vatory Sky Survey. The field is 6' x 6'. North is up, 
and east is to the left. The large and small circles are 
the WFC and NFI error circles, respectively. The opti- 
cal counterpart is indicated by tick marks. Its position is 
a 2 ooo.o = 22 h 39 m 20 s .90, Wo = +61°16'26".8, 



of emission from an unseen compact companion such 
as a neutron star. On this basis we identify the Be 
star with SAX J2239. 3+6116. This star is listed in the 
USNO-A2.0 catalog (Monet et al. 1996, Monet 1998), at 
a position a 20 oo.o = 22 h 39 m 20 s .90, <5 2 ooo.o = +61°16'26".8 
(uncertainty 0".3), which is 0.'3 from the NFI centroid 
and well within the error circle. The USNO-A2.0 catalog 
gives approximate magnitudes of B = 16.2 and R = 14.1. 

The summed spectrum from 3000 s of exposure is 
shown in Fig. [l^. We measure approximate magnitudes 
B = 16.5, V = 15.1, and R = 14.1 from the spectrum. The 
estimated uncertainty in these numbers is 0.2 mag. The 
magnitudes are consistent with the USNO magnitudes. 
The Ha emission line has an equivalent width (EW) of 
6.7 A, a flux of 3.4 x 10~ 14 ergs cm" 2 s -1 , and a full- 
width at half maximum (FWHM) of 400 km s -1 after 
correcting for the instrumental resolution. Emission lines 
of He I A5876, He I A6678, and He I A7064 appear to be 
double-peaked, with the peaks separated by 300 km s _1 . 
Double-peaked emission lines are commonly seen in Be 
star spectra, where they are attributed to a circumstellar 
disk. 
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Fig. 10. Optical spectrum of SAX J2239. 3+6116 obtained on the KPNO 2.1m telescope on 1999 December 2 and 3. 
In the lower panel, the flux scale refers to the upper trace of the spectrum. The lower trace is the same spectrum 
divided by a factor of 6. 



The reddening to the star may be estimated from many 
diffuse interstellar band (DIBs) that are apparent in the 
spectrum (Fig. [n]). We measure equivalent widths of EW 
= 3.4 A, 3.2 A, and 2.3 A for A4430, Na I D, and the 
AA5780,5797 blend, respectively. From the calibrations of 
Herbig (1975) and Tug & Schmidt-Kaler (1981) we esti- 
mate E(B — V) = 1.4 ± 0.2. This is less than the value for 
the total Galactic extinction as derived from the hydrogen 
column density measurements through HI maps (Dickey 
& Lockman 1995). The integrated H I column density is 
A^hi = 1.0 x 10 22 cm -2 in this direction. From the standard 
conversions Nn/Ay = 1.79 x 10 21 cm -2 mag^ 1 (Predehl 
& Schmitt 1995) and Ay = 3A3E(B - V), one obtains 
E(B — V) ~ 1.8. Also, the reddening is less than deter- 
mined from IRAS 100 /^m dust maps (Schlegel et al. 1998) 
which give E(B - V) ~ 2.0. 

Assuming that Ay — 3.13 E(B — V), the dereddened 
magnitudes from the spectrum become B = 10.7, V = 
10.7 and R = 10.9. Since these are somewhat redder than 



the colors of an early B star, we suspect that there may be 
some additional circumstellar extinction, and that dered- 
dened V ~ 9.7 is a realistic estimate. Also, if we assume 
an average Ay = 1 mag kpc" 1 , then the distance to 
SAX J2239.3+6116 can be estimated as 4.4 kpc. The ab- 
solute visual magnitude would then be -3.5, as expected 
for a star in the range BO V to B2 III. 

In a 16 square-degree area around the NFI position, 
there are 3 cataloged Be stars brighter than V = 10. None 
of these are coincident with the NFI error box. 

7. Association with other sources of high-energy 
emission 

In Fig. |ll], the WFC- and MECS-determined positions 
of SAX J2239. 3+6116 are shown, as well as those of three 
other possibly related X-ray sources. It appears likely that 
SAX J2239. 3+6116 is the same source as the one seen in 
the early seventies with Uhuru (4U 2238+60, Forman et 
al. 1978), OSO-7 (1M 2233+595, Markert et al. 1979), and 
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Fig. 11. Map of relevant high-energy sources. The large 
error circle to the right is that of 3EG J2227+6122, the two 
straight lines refer to the fast transient AT 2238+584, the 
dashed quadrilateral to 3A 2237+608, the solid quadrilat- 
eral to 4U 2238+60, the two small circles to the WFC 
error box (large) and NFI error box (small) of SAX 
J2239.3+6116. 



Ariel-V (3A 2237+608, Warwick et al. 1981). 4U 2238+60 
was detected at a level of (7 ± 1) x 10~ n erg cm~ 2 s _1 
(2-10 keV, statistical error only) which is about 4 times 
lower than the peak of the March 1997 outburst. This is 
not really surprising because neither measurement encom- 
passed the whole outburst. The Ariel source (Warwick et 
al. 1981) is classified as an irregular variable source during 
the 1974-1980 operational period with a measured peak 
flux of 2 x 10~ 10 erg cm~ 2 s _1 (2-10 keV) which is also of 
the same order of magnitude as SAX J2239.3+6116 and 
4U 2238+60. 

There are two other interesting sources which have er- 
ror boxes that are close to that of SAX J2239. 3+6116 but 
are formally inconsistent. The first one is the fast tran- 
sient AT 2238+584 which was seen in July 1975 with the 
sky-survey instrument on Ariel V (Pye & McHardy 1983). 
The outburst was not covered completely. The maximum 
measured flux was (1.6 ± 0.2) x 10~ 9 erg cm~ 2 s _1 (2-10 
keV) which is 5 times as bright as for SAX J2239.3+6116. 
The c-folding decay time of AT 2238+584 was 0.78 d. 
The other inconsistent source is 3EG J2227+6122 which 
is one of the unidentified EGRET sources (Hartman et 
al. 1999). We mention this because SAX J2239.3+6116 is 
the bright X-ray source closest to 3EG J2227+6122 and 
because SAX J2239. 3+6116 has quite a hard spectrum in 
the X-ray range. Nevertheless, the distance between the 
centroids of both sources is 87' while the 95% confidence 



error region of 3EG J2227+6122 is three times smaller 
(28'). 



8. Discussion 

SAX J2239.3+6116 is an X-ray transient which often re- 
curs with a periodicity of 262 d. Because of the Be-star 
nature of the likely optical counterpart, the periodicity 
may be identified with the orbital period of the binary. 
In such a Be X-ray binary system, the transient nature 
of the X-rays is thought to arise from a combination of 
episodic mass loss by the Be star, and an eccentric binary 
orbit. The compact object orbiting the Be star will show 
enhanced X-ray emission near periastron where it accretes 
more matter due to an enhanced wind density (e.g., review 
by White et al. 1995). The 2-10 keV flux history shows 
fluctuations of a factor of 10 3 . Possibly the true dynamic 
range is larger. The likely counterpart at an estimated 
distance of 4.4 kpc implies that the 2-10 keV luminosity 
ranged between lxlO 33 erg s^ 1 and 7 x 10 35 erg s^ 1 . It 
seems plausible that SAX J2239. 3+6116 was more or less 
active over the last decades. The peak flux is so low that it 
can easily be missed. The spectral type of the counterpart, 
the resulting luminosity and the transient behavior is con- 
sistent with the identification of SAX J2239. 3+6116 as a 
Be X-ray binary. This is also confirmed by the value for 
the X-ray to optical color index £ = B + 2.51og(Fx(/xJy)) 
(as defined by Van Paradijs & McClintock 1995) where 
B is the dereddened B magnitude and Fx the 2-10 keV 
X-ray flux. For B = 10.7 and F x = 0.03 /xJy (during 
quiescence) £ = 6.9. If B would be 9.7 (to correct for ad- 
ditional circumstellar extinction) then £ = 5.9. These are 
perfectly normal values for Be X-ray binaries in quiescence 
and three times lower than for low-mass X-ray binaries. 

Usually Be X-ray binaries contain an X-ray pulsar 
(e.g., review by White et al. 1995). SAX J2239.3+6116 
does not appear to contain one with a period shorter than 
10 3 sec, with a fairly sensitive upper limit of 10% on the 
amplitude. This may be explained by the fact that the 
pulse period is longer than 10 3 sec. Corbet (1984) discov- 
ered a strong correlation between orbital and pulse period 
for Be X-ray binaries. If SAX J2239. 3+6116 adheres to 
this correlation, the orbital period of 262 d implies a pulse 
period near to 10 3 sec. More dedicated sensitive X-ray ob- 
servations during future outbursts are obviously needed to 
study this system in more detail. 
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